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p Abstract 

In the MSSM, cosmological scalar field condensates formed along flat direc- 

- i— i i 

tions of the scalar potential (Affleck-Dine condensates) are typically unstable 
with respect to formation of Q-balls, a type of non-topological soliton. We con- 
sider the dynamical evolution of the Affleck-Dine condensate in the MSSM. We 
discuss the creation and linear growth, in F- and D-term inflation models, of the 
quantum seed perturbations which in the non-linear regime catalyse the collapse 
of the condensate to non-topological soliton lumps. We study numerically the 
evolution of the collapsing condensate lumps and show that the solitons initially 
formed are not in general Q-balls, but Q-axitons, a pseudo-breather which can 
have very different properties from Q-balls of the same charge. We calculate the 
energy and charge radiated from a spherically symmetric condensate lump as it 
evolves into a Q-axiton. We also discuss the implications for baryogenesis and 
dark matter. 
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1 Introduction 



Affleck-Dine (AD) baryogenesis [0] is a natural candidate for the origin of the baryon 
asymmetry |2| in the context of the MSSM and its extensions ||, with potentially 
important consequences for cosmology such as observable isocurvature perturbations 
non-thermal relic neutralinos || and number densities of baryons and dark mat- 
ter particles which are similar || [7|. The original picture of AD baryogenesis was 
of a homogeneous squark condensate formed along an F- and D-flat direction of the 
MSSM scalar potential, in which an asymmetry is induced via non-renormalizable soft 
SUSY-breaking terms. This subsequently thermalizes or decays, leaving the observed 
baryon asymmetry. More recently it was shown that the homogeneous AD condensate 
is typically unstable with respect to spatial perturbations || || . This is because the po- 
tential along flat directions of the MSSM scalar potential usually deviates from a pure 
(f) 2 potential, resulting in a negative pressure whenever the potential is flatter than </> 2 . 
Such deviations can occur because of A-terms in the potential || , because of a sharp 
change in the scalar mass term at large field values (gauge-mediated SUSY breaking 
models [R |HJ), or because of radiative corrections from the gauge sector (conventional 
gravity- mediated SUSY breaking models 0, 0]). It is the latter on which we will focus 
here. Although the main interest in unstable scalar field condensates along MSSM flat 
directions is perhaps AD baryogenesis, we emphasize that the phenomenon can arise 
for any coherently oscillating scalar field along a flat direction of the MSSM scalar 
potential with gauge interactions. This is because the gauge corrections will cause the 
mass squared term of the flat direction to decrease with increasing scale ||. In the 
present view of post-inflationary SUSY cosmology, the formation of such coherently 
oscillating scalar field condensates is a common and natural occurance, as a result of 
order H corrections to the SUSY breaking terms in the scalar potential |T1, I2|. In 
the following we will generically refer to the complex scalar field as the Affleck-Dine 
(AD) scalar. 

Spatial "seed" perturbations in the AD condensate, arising from quantum fluctua- 
tions of the AD scalar during inflation, will grow and eventually go non-linear, forming 
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condensate lumps || [/J. However, up until now there has been no clear picture of what 
happens to the lumps once they go non-linear. It has been assumed that they even- 
tually reach their lowest energy state for a given global charge, namely a Q-ball 
One aim of the present paper is to clarify the non-linear evolution of these condensate 
lumps. We will see that the condensate lumps generally do not initially form Q-balls, 
but rather a higher energy state which we refer to as a Q-axiton. These are somewhat 
similar to the original axitons JHJ (a form of pseudo-breather soliton \TE\) which form 
in the case of a real axion field, but in our case they can carry a charge by virtue of the 
complex nature of the AD scalar. Only in the case where the original AD condensate 
is carrying a near maximal charge are the properties of the Q-axitons close to those of 
Q-balls of a similar charge. 

For some of the important applications of Q-balls to cosmology, the ratio of the 
charge trapped in baryonic Q-balls (B-balls) to the total charge, fs, plays an important 
role |5|, [7], [16|] . A second goal of this paper is to estimate fs for various examples. 

The paper is organized as follows. In section 2 we discuss the cosmology of flat 
directions in SUSY models. In section 3 we discuss the quantum seed fluctuations. In 
section 4 we discuss the linear evolution of the spatial perturbations of the AD field. 
In section 5 we study numerically the non-linear evolution of a spherically symmetric 
condensate lump and estimate js for various examples. In section 6 we discuss the 
implications for AD baryogenesis and late-decaying Q-ball cosmology. In section 7 we 
summarize our conclusions. 



2 MSSM Flat Directions in Cosmology 

An F- and D-flat direction of the MSSM with gravity-mediated SUSY breaking has a 
scalar potential of the form || [7| 

where m is the conventional gravity-mediated soft SUSY breaking scalar mass term 
(m pa 100 GeV), d is the dimension of the non-renormalizable term in the superpoten- 
tial which lifts the flat direction, cH 2 gives the order if 2 correction to the scalar mass 
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(with c positive and typically of the order of one for AD scalars ||12|| ) and we assume 
that the natural scale of the non-renormalizable terms is M*, where M* = Mpi/8n 
is the supergravity mass scale [[|. The A-term also receives order H corrections, 
A x = Ax o + a x H, where A x is the gravity-mediated soft SUSY breaking term and ax 
depends on the nature of the inflation model; for F-term inflation \a\\ is typically of 



the order of one Jl2| whilst for minimal D-term inflation models it is zero [|17j . 



The logarithmic correction to the scalar mass term, which occurs along flat direc- 
tions with Yukawa and gauge interactions, is crucial for the growth of perturbations of 
the AD field and the formation of Q-balls |6|, [/J. This growth occurs if K < 0, which 
is usually the case for AD scalars with gauge interactions, since K is dominated by 
gaugino corrections 0. For flat directions with squark fields, K « —(0.1 — 0.01) || [7|. 
The actual value depends on the nature of the flat direction; in practice we will mostly 
be concerned with the case of baryonic Q-balls made of squarks. 

If cH 2 is positive at the end of inflation, the AD scalar will be at a non-zero 
minimum of its potential. Once H < m, the AD scalar will begin to coherently oscillate 
around its minimum at zero and the A-term will induce a global charge asymmetry 



in the condensate [12], jig] . The resulting asymmetry in the Universe at present is 



then proportional to the reheating temperature after inflation (Tr), 

2irn B T R 

VB = JPMk> (2) 

where tjb is the present baryon to entropy ratio and n# and H are the baryon charge 
density and the expansion rate when the asymmetry is formed (H « m) JF|. We see 
that, because of the dependence on Tr, the baryon asymmetry in AD baryogenesis 
cannot be determined independently of the inflation model; in practice, we can use 
the present baryon asymmetry to fix the reheating temperature and so determine 
the background cosmology. In general, the Universe will be matter dominated by 
coherently oscillating inflatons when the B asymmetry forms at H m m, since the 



thermal gravitino upper bound on Tr, Tr ^ 10 9 GeV |[L8||, implies that the value of H 



when the inflaton matter domination period ends at Tr is less than 1 GeV. 
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3 Seed Perturbations 



3.1 Growth of Perturbations 

Just like for the inflaton, there will be a spectrum of spatial perturbations induced in 
the AD scalar field by quantum fluctuations during inflation ||. For fields of mass 
m "Z H, the scalar fields have fluctuations on leaving the horizon of magnitude 

In general, the equations of motion for perturbations about the minimum of the po- 
tential have the form, 

5$ + 3H5<p - V 2 5(f) = -kH 2 5(j) , (4) 

where k is determined by the parameters c, \a\\ and d of the AD potential. For 
perturbations larger than the horizon, the spatial derivative can be ignored, since the 
perturbation can be treated as homogeneous on sub-horizon scales. During inflation 
(with H = Hi = constant), the solution of Eq.(4) for the amplitude of larger than 
horizon perturbations is 



oca- Re ^ ; ( 7=~ 



3 - V9 - 4fcl . (5) 



k > 9/4 gives a damped and oscillating solution for 8<p, for which the suppression of 
the amplitude of the perturbation by the expansion of the Universe is maximal, whilst 
k < 9/4 gives a purely damped solution. In general, Re(a) < 3/2. During the matter 
dominated period following inflation (with H oc a _3//2 ), the solution for the amplitude 
of 5(f) is 

1 
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(6) 



k > 9/16 gives a damped and oscillating solution for 5<p, with maximal amplitude 
suppression, whilst k < 9/16 gives a purely damped solution. In general, Re(rj) < 
3/4. Thus when the perturbation re-enters the horizon during matter domination, the 
amplitude is given by 

n X Re( v ) „-ite(<r)AjVe(A) Mr 



where AiV e (A) is the number of e-foldings before the end of inflation at which a per- 
turbation of scale A exits the horizon, a e is the scale factor at the end of inflation 
and a>hor is the scale factor when the perturbation re-enters the horizon. Once inside 
the horizon, the spatial derivatives become important. For ^ > H 2 and H > m, the 
spatial derivative dominates the potential term and <50 oc a -1 . Thus, writing in terms 
of the scale A at H during matter domination, the perturbation spectrum once the 
peturbation is within the horizon is given by 

5m W {H2 ^\ )X ; X = 1 " l Re {V) ~ l Re {a) ■ (8) 

For the case with the largest suppression of the amplitude (k > 9/4), Re (a) = 
2Re(r)) = 3/2, so that x = and the resulting spectrum is very simple^] 

WA) « ■ (9) 

An important feature of the spectrum in what follows is that the perturbation is 
usually larger for smaller values of A. 

3.2 Perturbations of the Affleck-Dine Field 

The actual value of k is determined by the form of the AD scalar potential. Let 
$ = (0i + «02)/v / 2 and let k\ and k% be the values of k for <50i and 502- in general, 
choosing the phase of $ such that XAx is real and negative, the minimum of the 
potential is at 02 = and 0i = 0i m , with 



c \ 2(d-2) Tr i \ 2 (d — 1) , , 

#3=3 ;a = V , 10 



where for simplicity we are considering the A-term along the 0i direction to be neg- 
ligible when minimizing; in general it results in a correction of order one to the 0i 
minimum. During inflation, H is constant and, on perturbing about the (0i m ,O) 
minimum, we obtain equations of the form of Eq. (4) with 

fci = c(2d-4) (11) 



1 This differs from the spectrum given in ||, which assumed a constant mass term throughout for 
the perturbation. 



for 5(f) i and 

k 2 = |a A |(d-l) 1/2 c 1/2 (12) 

for 5<p2- 

During the matter dominated period the minimum of the AD potential is time- 
dependent. As a result, when we perturb about (f>i m by A0i, we obtain 

Afa + 3#A0i = - ( d ~f H 2 ^ m - c(2d - 4)# 2 A0i . (13) 
4 [a — 2) 

Thus 0i will not oscillate about <f>i m , but rather about a larger value given by 

*® = { 1 + (14) 

Perturbing (0i, 2 ) about the (0i(£), 0) effective minimum then gives equations of the 
form of Eq. (4) with 

/«, ,n 9 (d- 3) (2d -3) 

fci = c(2d-4) + -^ ^ — ^ ^ 15 

V ; 4 (d-2) 2 V ; 

and 

Ojd-^ 1/2 /2 / 9(d-3M 

4(d-2) 2 I K ; V 8c (d-2) 2 / 1 ; 

If oa = 0, as in D-term inflation models, then the solution for the evolution of <50 2 is 
simply 5<p2 oc <f>i(t). This is easily understood, as in this limit there is no potential 
along the 9 direction of the AD field and so 59 = 50 2 /0i is constant. 

From this we see that, in general, the evolution of the 0i and 2 perturbations at 
H > m can be different, depending on c, d and a\. However, for c and \a\\ ~ 1, as 
expected in F-term inflation models, it is likely that fcj > 9/4 and so the amplitude 
of 5(f>i and £0 2 will be equally (i.e. maximally) suppressed throughout. In this case 
501 rs 50 2 at if m. We will focus on this case in the following. For \a\\ = 0, as 
expected in D-term inflation models, 50 2 oc 0i and so 59 is constant. In this case 50 2 
can be much larger than <50i at H ps m. 

To discuss the evolution of the perturbations at H < m, we need to know the 
initial spectrum of spatial perturbations of the AD field when the condensate and 
baryon asymmetry forms at H ps m. The directions 0i and 2 refer to the real and 
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imaginary directions as determined by the phase of the A-term. However, the phase of 
the A-term changes as H becomes smaller than to and A\ comes to dominate a x H. 
Therefore, 

(i) For the \a\\ ~ 1 case, corresponding to F-term inflation, the effect is to rotate 0i 
and 02 such that, if we consider the relative phase of A\ and a\H to be of the order 
of one (corresponding to the case of a condensate with amplitude 0i ~ 02 at H < to), 
then regardless of whether 50i or 50 2 dominates at H > to, 50i will be approximately 
equal to 502 once H < m. 

(ii) If |oa] ~ 1 but the relative phase is much smaller than 1 (corresponding to the 
case of a condensate with amplitude 0i ^> 2 at H < to), then it is possible for 502/02 
to be much larger than 50i/0i at if < to, for example if 50i ~ 50 2 &t H > m and 
0i ^ 02 at if ^ m. 

(iii) For the case with cza = 0, corresponding to D-term inflation, the initial phase of 
the AD field will initially be of order one relative to that of A Xo , so that, as with (i), 
50i ~ 502 and (pi 2 once H < m. 

Thus F-term inflation models can produce condensates which have different values 
for 50j/0j {i = 1,2) at H m, whereas minimal D-term inflation models always 
produce roughly equal values for 50i/0j. 

4 Linear Evolution 

Because of the logarithmic correction to the AD scalar potential, there will be an 
attractive force between the condensate scalars which causes the spatial perturbations 
to grow. The inflationary perturbations are known to have a small amplitude so that 
the initial growth will take place in the linear regime. Condensate fragmentation 
and collapse to soliton lumps happens in the non-linear regime, for which the linear 
evolution provides calculable initial conditions. Therefore we will first consider in 
detail the linear evolution of the perturbations. The homogeneous AD condensate can 
be characterized by the charge asymmetry it carries. The condensate is generally a 
sum of two real oscillating scalar fields. The maximum possible charge corresponding 
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to a given maximum amplitude of the AD scalar occurs when the combined oscillation 
corresponds to a circle in the (</>i,</>2) plane. We will refer to this as a maximally 
charged (MAX) condensate. This is roughly expected to occur in the case of F-term 
inflation models with an order one CP violating phase and in the case of minimal 
D-term inflation models. The condensate can also carry a less than maximal charge, 
as in the case of F-term inflation models with a small CP violating phase, in which 
case it will describe an ellipse in the (<pi,<p2) plane. 

4.1 MAX Condensate 

A solution for the linear evolution of spatial perturbations of the condensate has been 
given previously for the case of a MAX condensate |7], §], which we refer to as the 
Kusenko-Shaposhnikov (KS) solution ||. It assumes a solution of the form <fr = <j)(t) + 
5<fi(x,t) and 9 = 6{t) + 59(x,t), where the homogeneous MAX condensate is described 
by 

$ = #& e m , (17) 

with cf)(t) = (a /a) 3 / 2 (j) (a is the scale factor) and 9(t) 2 « m 2 , up to corrections of 
order Km 2 . (As it is true for most models, we will assume that \K\ is small compared 
with 1.) The KS solution also requires that S(f)(x,t) and 56(x,t) initially satisfies 

-(f),- <*> 

As discussed in the previous section, for a MAX condesate we expect 6(pi/4>i ~ 502/02, 
so that this condition is satisfied. The solution of the linear perturbation equations 
then has the form |7| 

-•".^f/.wV",^ (19) 



and 



59 « SB, ex P | / dt (^^) 1/2 )e ikx • (20) 

These apply if |k 2 /a 2 | ^ \2Km 2 \, H 2 is small compared with m 2 and \K\ <C 1. If 
the first condition is not satisfied, the gradient energy of the perturbations produces 
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a positive pressure larger than the negative pressure due to the attractive force from 
the logarithmic term, preventing the growth of the perturbations. 

For the case of a matter dominated Universe, the exponential growth factor is 

, flk 2 \K\m 2 \ l/2 2 /|K|k 2 \ 1/2 , . 

dt k-7^T =TF > (21) 



,2 a 2 9(t) 2 J H V 2 a 2 

where we take the scale factor when the AD oscillations begin to be equal to 1. The 
largest growth factor will correspond to the largest value of k 2 for which growth can 
occur, k 2 /a 2 ~ 2\K\m 2 . Thus the value of H at which the first perturbation goes 



non-linear is 



with 



H, « «=! , (22) 
a{A) 



Q .(A) = -log (*M*>) , (23) 

where (f) is the value of <fi when the condensate oscillations begin at H ~ m. A typical 
value of a (A) for the spectrum Eq.(9) (e.g. for d — 6) is a (A) ~ 30. The initial 
non-linear region has a radius \ at Hi given by 

Thus for a MAX condensate, when the condensate is just going non-linear, the ampli- 
tude of the scalar field is given by 

0(x, t) = <p(t) + <ty(x, t) w 0(t)(l + cos(k.x)) , (25) 

where |k| is given by the wavelength which first goes non-linear, Aj. This will provide 
the initial condition for the numerical study of the non-linear regime. 

4.2 Non-MAX Condensate 

For the case of a non-MAX condensate, we cannot directly use the KS solution. How- 
ever, it is likely that the initial radius and the time at which the spatial perturbations 
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initially go non-linear will roughly be the same as for the MAX condensate. We will as- 
sume that 5(pi ~ 5<p2, as is true if both scalar perturbations are maximally suppressed 
when H > m. The equations of motion for 0j (i=l,2) are given by 

0. + _ V 2 0* = -m 2 (l + K)<j) t - Km 2 Uog (^^j • (26) 

For 0! ^> 2 , the equation for 0i will be similar to the case of the MAX condensate, 
except the 2 + (ft 2 , will be oscillating in time rather than constant. Thus the equation 
for the growth of perturbations in X will be similar to the MAX condensate and the 
0i perturbations will go non-linear (<50i/0i > 1) roughly as in the case of the MAX 
condensate. The equation for perturbation in 2 will, however, be different, since the 
log term is dominated by <p\. The condition for the 02 equation to go non-linear is 
then that 502 ^ 0i. We expect the rate of growth of (502 to be no larger than for 5 (pi, 
since the perturbations grow due to the attractive interaction between the scalars due 
to the log term, and the interaction of 2 with log((pf) is the same as that of 4>i- So 
if 501 ~ 502 initially then we expect (502 ^ (50i throughout. Therefore non-linearity 
will occur only once <50i goes non-linear, at which point the condensate will begin to 
fragment to condensate lumps. In general, the charge density of the initial non-linear 
lumps will be essentially the same as that of the original homogeneous condensate. 



5 Non-linear Evolution and fs 
5.1 Initial Lump 

In order to study the non-linear evolution, we will consider the evolution of a single 
condensate lump. We can neglect the expansion of the Universe, since at the time when 
the spatial perturbations go non- linear, the expansion rate, Hi, is small compared with 
the dynamical mass scale in the equations of motion governing the growth of the lumps, 
\K\ 1 / 2 m, so that the 3i70 terms in the equations of motion are negligible compared 
with the other terms. From the discussion of Sect. 4, when the spatial perturbation is 
just going non-linear (5(f) 0), the AD field for the MAX condensate can be written 
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as 

Ae imt 

$(x,t)«-^(l + cos(k.x)), (27) 

where A/ y/2 is the amplitude of the coherent oscillations when the field first goes non- 
linear. In terms of $ = (0i + i(f>2)/y/2, the initial non-linear perturbation is described 
by 

0i(x,i) = Acos(mt)(l + cos(k.x)) (28) 

2 (x, t) = Asin(mt)(l + cos(k.x)) . (29) 

This can be thought of as a 'lattice' of adjacent condensate lumps. In order to gain 
some insight into the evolution of these lumps, we will consider a single spherically 
symmetric lump. Such lumps are described in general by 

<f>i(r,t) = Acos(mt)(\ + cos(irr/2r )) (30) 

< t>2(j'it) = Bsin(mt)(l + cos(nr /2r Q )) , (31) 

for r < 2ro and by 0^2 = otherwise. The initial radius of the lump is 2tq, where 
r = n /(y^lK^^m). The MAX condensate lump corresponds to A = B, whilst the 
non-MAX lump has A > B. The corresponding energy and charge densities (with unit 
charge for the scalars) are 

p= |$| 2 + |V$| 2 + V(|$|) (32) 

and 

q = X 2 - X 2 , (33) 

where 



1$ 



2 X 



\/(|$|)=m 2 |$| 2 h+Klogl^) ) . (34) 



The total energy and charge in a fixed volume are given respectively by 

E = 4tt J drr 2 p , Q = 4vr J drr 2 q ~ AB, (35) 
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with, for the initial MAX condensate lump of Eqs.(30, 31), 

Q = Qmax = 132.2A 2 \K\~ 3/2 m- 2 . (36) 

The non-maximally charged lump with B < A then corresponds to Q/Q maiX = B/A < 
1. 

5.2 Numerical Solution 

The scalar field equations of motion are given by Eq.(26) with H = 0. We have solved 
these equations numerically for the case of the spherically symmetric lump. We choose 
A = 1 and set m = 100 GeV when needed. To avoid the singularity at \<p\ = 0, where 
the one-loop logarithmic correction alone no longer is valid, we have introduced a cut- 
off 5 by letting (<f)\ + (fil)/(fil — > {4>\ + 0i)/0o + $\ we nave verified that the value of the 
small cut-off does not significantly affect the solutions. With the initial lump radius 
2r we have considered a spatial sphere of radius 8r , at the boundary of which the 
outflowing waves are damped by hand. This is to prevent reflected waves bouncing 
back onto the lump, a situation which would be realistic only if the average lump 
distance were extremely small. 

The behaviour of the solutions depends on K, and to a greater extent on Q/Q max . 
In general, the condensate lump pulsates while charge is flowing out until the lump 
reaches a (quasi-)equilibrium pseudo-breather configuration^ (in which the lump pul- 
sates with only a small difference between the maximum and minimum field ampli- 
tudes), as seen in Fig. la. We refer to this state as a Q-axiton. For non-maximal 
condensates, it is very different from the lowest energy configuration for a given charge, 
namely a Q-ball. The Q-ball is essentially made entirely of charge, with E w mQ (ne- 
glecting the small binding energy per charge) || [/J. For the Q-axiton, in which the 
attractive force between the scalars is balanced by the gradient pressure of the scalar 
field, the energy per unit charge can be much larger than m; indeed, the Q-axiton 

exists even if Q = 0. Only for a more or less maximally charged Q-axiton are the 
2 We distinguish between spatial pulsations of the lump as it collapses and breathing, which is 
essentially the coherent oscillation of the AD field within the lump. 
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Figure 1: (a) The time evolution of the field amplitude at the origin, 0i(O); (b) a detail 
of the early time evolution, showing the pulsation cycles modulating the coherent 
oscillations (time is in units of 1/m); (c) the first pulsation of the condensate lump: 
a) initial lump; b) maximal lump; c) minimal lump; (d) the final equilibrium Q-axiton 
compared with the initial lump (r is in units of ro). In all cases K = —0.05 and 

Q Qmax- 

properties close to that of the corresponding Q-ball. We believe that the Q-axiton 
should eventually evolve to the lower energy Q-ball stateQ, although this is not appar- 
ent from our numerical results and we cannot rule out the possibility that the Q-axiton 
may be very long-lived. Compared with the natural time scale t = 1/m, the time taken 
to reach the Q-axiton state is long, of the order of 1600/m for the case K = —0.05 
and Q = Q ma , x , which is about five expansion timesQ (H~ l rs 300/m). This is shown 
in Fig. 1, where we display the oscillation of the field amplitude at the origin as a 
function of time, as well as the spatial profile of the whole lump as it pulsates. Fig. 

3 There are no absolutely stable spherically symmetric breather solutions in 3 + 1 flat space |p^| . 
4 Expansion can, however, be neglected in numerical solution for the Q-axiton itself, as noted 

previously. 
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Figure 2: (a) The time evolution (in units of 1/m) of the Q-axiton energy E and (b) 
charge Q for K = —0.05 and initial Q = Q max . 

la shows the rapid pulsations gradually settling towards the equilibrium by emitting 
scalar field waves. Fig lb shows the first four pulsations of the lump; the high fre- 
quency oscillations modulated by the pulsations are the coherent oscillations of the 
AD field. The lump keeps its shape during the whole time while pulsating and slowly 
evolving to the final configuration, as depicted in Figs, lc and Id, which show the 
maximum amplitudes at different stages of the pulsations. In the case with Q = Q ma x, 
4>i and (f>2 have almost identical time evolution, the only difference being that there is a 
phase difference of it/ 2 between the 4>i and 2 oscillations, as can easily be understood 
from the symmetry of the initial conditions and scalar field equations of motion. 

In Fig. 2 we show the time evolution of the charge and the energy of the whole 
configuration, integrated out to the distance 8r . At first charge and energy flows out 
of the volume, but the slow approach to equilibrium can also readily be seen, with 
the axiton lump starting with initial charge Q = 1.18 and stabilizing to Q ~ 1.09. 
Q and E are nearly proportional, as they should be for the case with Q m ax and <pi 
and 02 different only by a phase tt/2. (In fact, E/Q decreases from 106 GeV initially 
to 98 GeV, due to the increased binding energy of the Q-axiton state.) We see that 
for the case of the MAX condensate the ratio of E to Q is close to m = 100 GeV 
throughout, showing that the Q-axiton formed in this case is essentially a Q-ball. 

In Fig. 3 we show the field amplitude and lump oscillations for the non-maximally 
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Figure 3: (a) The time evolution of the field amplitude at the origin, 0i(O); (b) a detail 
of the time evolution, showing both 0i(O) and 02(0); (c) a further detail. In all the 
cases K = —0.05 and Q = O.lQmax while time is in units of \ jm 

charged case Q = 0.1Q max . Now the squared 0i an d 02 amplitudes do not follow 
a circle as in the maximally charged case, but rather a precessing ellipse (Fig. 3b). 
The qualitative behaviour is similar to the maximally charged case in that while E 
and Q at first escape the fiducial volume, at later times the configuration reaches a 
quasi-equilibrium state (Fig. 4). The initial fluctuation in Q is a numerical artifact; 
in general, our numerical solution for the evolution of the charge of the lump becomes 
less accurate for smaller Q/Q m ax- The ratio of E to Q of the Q-axiton in this case is 
approximately 600 GeV, much larger than for the corresponding Q-ball. 

5.3 Values of f B 

We have also calculated the ratio of the charge in the equilibrium Q-axiton state 
to the charge of the initial condensate lump. The results as a function of Q and K are 
given in Table 1, from which we see that decreases for smaller values of Q and for 
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Figure 4: (a) The time evolution (in units of 1/m) of the Q-axiton energy E and (b) 
charge Q for K = —0.05 and initial Q = O.lQmax- 

larger vales of \K\. 

Table 1. Values of f B 



K 


Q / Qmax 


fB 


-0.1 


1 


0.89 


-0.05 


1 


0.92 


-0.01 


1 


0.94 


-0.1 


0.1 


0.64 


-0.05 


0.1 


0.74 


-0.01 


0.1 


0.83 


-0.1 


0.01 


~ 0.3 


-0.05 


0.01 


~ 0.5 


-0.01 


0.01 


~ 0.6 



In all this we have considered a single spherically symmetric lump isolated in space, 
compared with the real case in which there is a lattice of closely spaced lumps which 
are in general not spherically symmetric. It is beyond the scope of the present paper 
to study this more realistic case, but we note that the rate at which the lumps radiate 
charge and energy is very slow. By the time the lumps have radiated their charge 
and reached the Q-axiton state, they will have been substantially pulled apart by 
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expansion, so that the lumps are eventually isolated. Therefore the results for for 
the spherically symmetric lump are likely to be a reasonable estimate of the values in 
the realistic case. Indeed, since we might expect non-spherically symmetric lumps to 
be less efficient in reaching their quasi-equilibrium configurations, the values of for 
the spherically symmetric lump may well be an underestimate of the realistic value. 
We will assume in the following that the Q-axitons eventually evolve to lower energy 
Q-balls of the same charge. To do this they must lose their excess energy, presumably 
by a very slow radiation of scalar field waves. However, the full evolution of the Q- 
axiton to the final Q-ball requires a much more ambitious numerical analysis than that 
attempted here. 



6 Consequences for Affleck-Dine Baryogenesis 

We have previously discussed the possibility that Q-balls can decay at a very low 
temperature, < 1 GeV. Such low decay temperatures naturally occur for the case 
of d = 6 AD baryogenesis, which has a reheating temperature typically around 1 GeV 
[0, Late decaying Q-balls ("B-ball baryogenesis" || [7J) can both protect the baryon 
asymmetry from the effects of L violating interactions (such as arise with Majorana 
neutrino masses) and allow for an understanding of why the number density of baryons 
and dark matter particles are within an order of magnitude of each other when the 
dark matter particles have masses of the order of raw- I n the case where the similarity 
of the number densities is explained by the late decay of Q-balls to lightest SUSY 
particles (LSPs) and baryons, the value of fs immediately gives us the mass of the 
LSP. The ratio of the number densities is given by J5], |?J 

n B 1 



(37) 



The observed range of ratios is [16 



(*L) =(1.5-7.3)^. (3* 

\n DM J nhs mw 



f^DM / obs m W 

To be consistent with present experimental limits on MSSM neutralinos, m x > 30 GeV 
PH[ ], we would require that ^ 0.6. From Table 1 we see that this rules out MSSM 
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neutralinos as the LSP for K in the range —0.1 to —0.01 and Q/Q ma x > 0.1. How- 
ever, for Q/Qmax = 0.01 it is probable that, for sufficiently large \K\, light MSSM 
neutralinos from Q-ball decay can be consistent with experimental bounds. Our nu- 
merical results become less accurate for smaller values of Q and larger \K\, but we 
expect that MSSM neutralinos up to 60 GeV can be consistent with experimental 
limits for K = —0.1 and Q/Qmax = 0.01. It is interesting to note that a recent dark 
matter search experiment has suggested the existence of a candidate of mass 59 GeV 



21 . It should be kept in mind that in the realistic non-spherically symmetric case fs 



could well be smaller, so allowing for larger neutralino masses. On the other hand, if 
fs ^ 0.6, then either we must allow the neutralinos to annihilate after Q-ball decay, 
so breaking the direct connection between the number densities of baryons and dark 
matter particles 0, [7j], or we must go beyond the MSSM. A light singlino in the next- 
to-minimal SUSY Standard Model (NMSSM) would be an interesting possibility for 



the LSP in this case [22 



This all assumes that the Q-axitons evolve into the corresponding Q-balls. Should 
the Q-axitons be very long-lived, their cosmology, for the case of a non-MAX conden- 
sate, could be substantially different from that of the Q-balls . 



7 Conclusions 

We have considered the origin and linear evolution, in the context of SUSY F- and D- 
term inflation models, of spatial perturbations of an Affleck-Dine condensate and the 
collapse of a spherically symmetric condensate lump of the type expected to arise from 
the fragmentation of the condensate. We find that for a typical range of parameters, 
30-90 % of the charge of the collapsing lump ends up in the final quasi-equilibrium 
state (fs ~ 0.3 — 0.9), with the remainder being radiated in the form of ripples of 
scalar field. The quasi-equilibrium state is, in general, not a Q-ball, but a higher energy 
pseudo-breather state, a Q-axiton. This physically approximates a Q-ball only for the 
case of an initial condensate which is nearly maximally charged. Assuming that the 
Q-axiton evolves into the corresponding Q-ball, for the case where dark matter and 
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baryons originate directly from late-decaying Q-balls, so explaining the similarity of 
the number densities of baryons and dark matter particles, fs ^ 0.6 allows for a range 
of neutralino masses consistent with experimental constraints on MSSM neutralinos. 
Such values can occur if the condensate is not maximally charged, fg > 0.6, on 
the other hand, requires rriLSP ~ 30 GeV, which rules out MSSM neutralinos coming 
directly from late-decaying Q-balls as dark matter, although an NMSSM singlino could 
be consistent with experimental constraints. The numerical results presented here are 
based on the evolution of a single, isolated, spherically symmetric condensate lump. 
To study Affleck-Dine condensate collapse for the realistic case, with many condensate 
lumps of different shapes and sizes, and to study the evolution of the Q-axitons into Q- 
balls, we would require a much more ambitious numerical analysis than that presented 
here. We hope to develop this in the future. 
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